Model investigation on the probability of QGP formation 
in relativistic heavy ion collisions * 
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The formation probability of quark-gluon plasma in relativistic heavy ion collisions for colliding 
nuclei of different sizes is investigated in the framework of a bond percolation model. The results 
show that nuclei with sizes smaller than that of Pb or Au produce QGP with probability less than 
unity even at very high collision energies. The dependence of the QGP-formation probability on 
different nuclear sizes and on various centralities of Au-Au collision are presented. 

PACS numbers: 25.75.Nq, 12.38.Mh, 64.60.ah 



I. INTRODUCTION 

It is believed that a medium with deconfmed quarks 
and gluons as constituents, referred to as quark-gluon 
plasma (QGP), could be created in high energy nucleus- 
nucleus collisions. Up to now, abundant experimental 
data have been obtained through the collisions of Si, O, 
Al, Cu, Au and Pb et al. at center-of-mass energies vary- 
ing from about 2 to 200 AGeV PUS, Hi, 0,0, II- The 
recent data analyses on RHIC Au-Au collisions at 200 
AGeV show strong evidences for liberated quark degree 
of freedom over nuclear volume d, 0, 0] . 

The widely accepted definition for QGP is: a (locally) 
thermally equilibrated state of matter in which quarks 
and gluons are deconhned from hadrons, so that color 
degrees of freedom become manifest over nuclear, rather 
than merely nucleonic, volumes (5j. From this defini- 
tion, formation of QGP requires quark deconfinement in 
a large volume at least larger than that of a nucleon. 
Therefore, proton-proton collision could not form QGP 
even at very high energy. Then it is natural to raise 
a question: what is the dependence of QGP-formation 
probability on the size of colliding nuclei? 

The above question could be understood qualitatively 
from the point of view of energy density deposited in 
the collision region. People believe that QGP could be 
created if the energy density is above about ten times 
that of the normal nuclear matter 0, [l^, HH. The 
energy deposition in collision region is due to multiple 
scattering between nucleons from two incident nuclei. 
In the Glauber model [12], which describes the multi- 
ple scattering of nucleus-nucleus collision, the probabil- 
ity of n multi-scattering in the collision of two nuclei 
at a given impact parameter is a binomial distribution. 
For two head-on nuclei with identical nucleon number A, 
the mean number of multi-scattering is proportional to 
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A 4 / 3 [l3|. For small nuclei the multi-scattering number 
n may be large enough for QGP formation only at the 
tail of n-distribution. So it could be inferred that if the 
colliding nuclei are not large enough, the probability for 
producing QGP in their collision will be much less than 
100%. 

In this paper we suggest to investigate QGP-formation 
probability in heavy ion collisions by using a bond per- 
colation model [H| . Application of percolation theory to 
quark deconfinement was first suggested b y G . B aym llal 
and further extended by H. Satz et al. p^.Tl^[Tl.ll9l.[2o|. 
In their work, they use site percolation model and dis- 
cuss the critical nucleon density of phase transition. In 
Ref. [HI the bond percolation model is applied to dis- 
cuss the cluster formation in analytic crossover between 
hadronic gas and QGP. We are now applying this model 
to study the probability of QGP formation. 

Firstly, let us have an intuitive look at the process that 
occurs in heavy ion collisions. When two nuclei collide 
with high velocity, due to Lorentz contraction the scale 
in the longitudinal direction is much less than that in the 
transverse plane, and the two nuclei can be described as 
two discs without thickness. During the collision, the 
nucleons in the two discs interact with each other and 
the potential barriers between them decrease as the in- 
crease of the center-of-mass energies. That is to say, the 
wave- function of nucleons will be distorted at high center- 
of-mass energies, and the infinitely high confinement- 
potential between neighboring nucleons might be reduced 
to a finite-height potential barrier, cf. the central sub- 
figure in Fig. 1(a). The higher the center-of-mass energy 
is, the more distorted is the nucleon wave function and 
the lower are the potential barriers. Because of quantum 
tunneling effect, quarks in nucleons are able to delocal- 
ize from single nucleon and the nearer the two nucleons 
the larger probability of derealization. Let us use S 
to denote the distance between the neighboring nucleons 
that have quark-delocalization, cf. the central sub-figure 
of Fig 1(a). At fixed energy, there exists a maximum 
distance S m!ix , such that quark derealization is possi- 
ble when S < S max , but is impossible when S > S max - 
As the increase of energy, the maximum derealization 
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FIG. 1: (Color online) A schematic plot for the bond forma- 
tion, (a) Nucleons connected by bonds form clusters. The 
central sub-figure shows how the infinite-high confinement 
potential between two neighboring nuclei reduces to a finite- 
height potential barrier. Quark tunneling through the barrier 
forms bond, shown as full line-segments in the main figure, 
(b) A cluster extending from one boundary to another is an 
infinite cluster. 



distance S max increases. 

If quark derealization happens between two nucleons, 
the two nucleons are connected by bond, shown as full 
line segments in Fig's. 1, to form cluster, inside which 
quarks are free to tunnel from one nucleon to the other 
and the nucleons turn to colored objects, referred to as 
cells. The size of a cluster is defined as the number of 
cells included in it. It can be seen from Fig. 1(a) that the 
clusters can be of various sizes. 

As the size of a cluster grows to the nuclear scale, i.e. 
extending from one boundary to the other, cf. the big 
cluster in Fig. 1(b), color degree of freedom are mani- 
fested over nuclear volume and QGP forms. 

The above picture could be realized by a two- 
dimensional bond percolation procedure. In the usual 
geometrical percolation model, the control parameter is 
the bond or site occupation probability p [2l[. When 
p equals a critical value, an infinite cluster appears and 
the system turns from disconnected phase to connected 
phase. In our bond percolation model, the control pa- 
rameter is the maximum derealization distance S max , 
which depends on the collision energy of the two nuclei. 
If two nucleons depart with distance less than S max , there 
could be a bond formed between them, representing the 
tunneling of quarks through potential barrier. The nucle- 
ons (cells) aggregate through bond-connection and form 
clusters. The appearance of infinite cluster changes the 
medium from color insulator to color conductor, or from 
hadron phase to quark-gluon phase. 



II. THE CONSTRUCTION OF THE BOND 
PERCOLATION MODEL 

In the bond percolation model, two head-on colliding 
nuclei with nucleon number A are simplified as two over- 



lapped discs of radius R = 1.2 A 1 / 3 and with 2 A nucleons 
randomly distributed inside the overlapped region. Each 
nucleon has a hard-core radius r e = 0.1 fm. A nucleon, 
or cell after cluster formation, departs from the center of 
the big disc father than R—r e is referred to as a boundary 
cell. 

The percolation procedure is as follows: 

(1) Randomly select a cell a as a mother cell. 

(2) Find the cells that can form bonds with the mother 
cell, which will be referred to as bond- candidate cells and 
are defined as those cells with \r — r a \ < S. It is assumed 
that each cell can be connected by 3 bonds at most since 
each nucleon has 3 constituent quarks. So we randomly 
select 3 cells from the bond-candidate cells to form 3 
bonds connected to the mother cell a. These are referred 
to as daughters. If the number of candidate cells is less 
than 3, then the number of daughters is equal to the 
candidate number. 

(3) Find the bond-candidate cells for the daughters 
of cell a. For every daughter find her bond-candidate 
cells from the remaining unbounded cells, and randomly 
select 2 bond-candidate cells to form bonds. The cells 
connected to daughters are called granddaughters. 

(4) Repeat the procedure to granddaughters and 
granddaughters' daughters we will get a cluster, which 
grows until no bond-candidate cell can be found anymore. 

(5) Then choose another cell (3 from the left unbounded 
cells as another mother cell, and repeat the procedure 
starting from step 2. 

In this way, every cell is assigned to a cluster. In ev- 
ery cluster, find the boundary cells if any, calculate the 
distance between every two boundary cells, and denote 
the maximum distance by d. A cluster with d > y/2R 
is called an infinite cluster. The probability Poo for the 
appearance of infinite cluster is defined as: 



where jY^ is the number of events with infinite cluster, 
,jV is the total number of events in the sample. In this 
model, Poo is the formation probability of QGP and will, 
therefore, be denoted in the following by Pqgp- 



III. THE DEPENDENCE OF THE 
PROBABILITY OF QGP FORMATION ON THE 
SYSTEM SIZE AND CENTRALITY 

The bond percolation simulation is done for nucleus- 
nucleus collisions of different nuclei. The variation of 
-Pqgp as a function of S max is shown in Fig. [2] for nuclei 
with different nucleon number A. It can be seen that 
for each kind of nucleus as the increase of S max , -Pqgp 
gradually increases from to a saturation value. This 
is typical for finite size percolation model, while for an 
infinite system .Pqgp (Smax) would be a step function and 
the point where Pqgp starts to be greater than would 
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FIG. 2: (Color online) The probability Pqgp of infinite cluster 
as a function of the maximum percolation distance S'max for 
nucleus- nucleus collisions of different nuclei calculated from 
bond percolation model. The statistical errors are within the 
symbols. The curves represent the fits to Eq. ([2]). 



FIG. 3: The dependence of maximum (saturation) QGP- 
formation probability -Pqgp on nuclear size (bottom scale) 
and on the centralities of Au-Au collisions at ^/snn = 200 
GeV (upper scale). 



be the threshold S c . In our case the system is of finite 
size, so we use a function p2| 

PQGp(SWax) = o[l + tanh[6(S* max - c)}], (2) 

to fit the shape of Pqgp (5 max ) in Fig. [2l where a, 6, c are 
parameters. It turns out that the fits are very good, cf. 
the curves in Fig. 

The inflexion point c of the fitting curve can be used 
as an evaluation of the threshold of S max : S c — c, and 
the saturation value of Pqgp is Pqgp = 2a. The re- 
sults obtained from the parameters of the fits are listed 
in Table [T] It can be seen from Table U that comparing 
with small nuclei, the bigger ones have higher saturation 
values PA GP and smaller S c . 



TABLE I: The saturation value Pqgp or Pqgp and the critical 
percolation distance S c for different size of nuclei. 





U 


Pb 


Au 


Sn 


Cu 


S 


Si 


O 


C 


A 


238 


207 


197 


119 


64 


32 


28 


16 


12 


Pqgp (%) 


100 


99.8 


99.8 


98.7 


91.4 


71.6 


66.6 


45.4 


35.6 


Sc (fm) 


0.67 


0.69 


0.70 


0.78 


0.90 


1.05 


1.08 


1.22 


1.29 



The Pqgp is the maximum (saturation) QGP- 
formation probability and the S c is the point where 
the QGP-formation probability is half of the saturation 
value. Both of them depend on the function Pqgp (■S'max) j 
where <!5 m ax is the maximum derealization distance. As 
the increase of center-of-mass energy ^snn, the dereal- 
ization will happen in an increasing range and »S* ma x in- 
creases. From Table HI we see that larger nuclei have 
smaller S c , which means that the energy threshold to 
form QGP for larger nuclei is lower than that of the 
smaller ones. The maximum QGP-formation probability 
Pq GP for smaller nuclei are lower than those of the bigger 
ones. For small nuclei the QGP-formation probability is 
less than 100% even at very large S max , or equivalently 



at very high center-of-mass energy ^/snn, £•&. the QGP- 
formation probability gets saturated. It can be seen from 
Table |T] that the saturation values of Pqgp for U-U, Pb- 
Pb and Au-Au collisions are about 100%, while those of 
smaller nuclei do not reach 100%. For example, Cu-Cu 
collisions only have about 91% events to form QGP at 
very large S max (very high ^/snn)- From the variation 
of Pqgp on nuclear size A shown in Fig. [3J we see that 
the saturation value of Pqgp decreases quickly when the 
nuclear size is less than that of copper. 

We further transform the dependence of maximum 
(saturation) QGP formation probability on nuclear size 
to that on the centrality of Au-Au collisions at ^/snn = 
200 GeV by using the function between centrality and 
number of participants calculated from the Glauber 
model [231 ] . The dependence of QGP-formation proba- 
bility on centrality in Au-Au collisions is also shown in 
Fig. [31 We see that for Au-Au collision at ^/snn = 200 
GeV, almost all the events with — 40% centralities could 
produce QGP, while in the most peripheral collisions (80- 
100% centrality) there are only about 40% events that 
can form QGP. 



IV. CONCLUSION AND DISCUSSION 

It is argued that since the number n of multi-scattering 
in high energy nucleus-nucleus collisions has a binomial 
distribution with the average value of n proportional to 
A 4 / 3 of the colliding nuclei, when the colliding nuclei are 
small the energy density of the medium produced in each 
collision may not be high enough to form quark-gluon 
plasma. Therefore, the probability Pqgp for QGP for- 
mation in high energy heavy ion collisions may not be 
100% in all cases. 

The bond percolation describes the quark derealiza- 
tion and transition from hadronic matter to quark-gluon 
plasma reasonably. Assuming that the production prob- 
ability of infinite clusters in the percolation model is the 
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formation probability of quark-gluon plasma, we found 
from the percolation simulation that the maximum (sat- 
uration) QGP-formation probability depends on the nu- 
clear size and collision centrality. Due to this depen- 
dence, when comparing experimental results of not very 



large nuclei and/or not very high centralitics with the- 
oretical predictions, an efficiency correction would be 
needed, since the data are a mixing of events with QGP 
and those without QGP. 

The authors thank Dr. Liu Hui for helpful discussions. 
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